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The fate of heavy quarkonia states has been investigated in QCD at high temperature when 
the plasma has a small momentum space anisotropy within a quasi-particle model. A real time 
static potential has been obtained from a Hard thermal loop expression for the gluon self-energy, by 
employing the quasi-parton equilibrium distribution functions extracted from hot QCD equations 
of state of 0{g'^) and 0{g'' ln(l/g)). Employing the potential, quarkonia melting has been studied. 
It is found that the equations of state and the momentum anisotropy cause significant modifications 
in the screening properties of the plasma. 
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I. INTRODUCTION 

The dissociation of heavy quarkonia states due to the 
presence of screening of static color fields in hot QCD 
medium has long been proposed as a signature of a de- 
confined medium, and QGP formatioii[l| . Since than, 
this has been a area of active research and undergone 
several refinements[l, 0, However, a precise under- 
standing of this phenomenon is still a subject of debate. 
There are disagreements between various approaches due 
to several ambiguities . Attempts are being made on 
to fix them[3]. A precise definition of the dissociation 
temperature is still elusive and is a matter of intense the- 
oretical and phcnomenological investigations -from the 
perspective of lattice spectral function studies 0, [TTI - [l3 | 
or potential inspired models [l3, [isl - fTTj or effective field 
theories [isj. Experimentally, quarkonia suppression has 
been seen in RHIC experiments at BNLQ and hopefully 
the understanding will be refined at LHC heavy ion pro- 
gram in CERN in the near future. 

The heavy quark-antiquarks, such as cc are bound to- 
gether by almost static gluons [l9l - [2l| . The fate of qq 
bound states, in a QGP medium, can be explored by 
determining the behavior of gluon self energy at high 
temperature. Recall that the real part leads to Dc- 
bye screening on which the traditional understanding of 
quarkonia melting in QGP medium is based. The imag- 
inary part of the potential implies that thermal effects 
generate a finite width for the quarkonium peak in the 
dilepton production rate. Thev|2l| further showed that 
the physics related to the finite width originates from 
the Landau-damping of low-frequency gauge fields, and 
could be studied non-pcrturbatively by making use of the 
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classical approximation. This leads to a new proposal to 
explain the quarkonia melting in QGP medium j2ll. [2^. 
According to which the quarkonia bound states dissoci- 
ate effectively at lower temperatures as compared to that 
obtained from Debye screening picture. In this picture, 
the quakonium state will disappear whenever the ther- 
mal width exceeds the binding energy. The temperature 
lies in the range ff^M < T < gM, where M is the heavy 
quark pole massyjj. The estimate has been refined em- 
ploying the understanding from the effective theories in 
[23l - l25l |. These studies suggest that the quarkonia disso- 
ciation in the thermal medium is not dominated by the 
Debye screening rather by the above mentioned mecha- 
nism. Further, modifications to these estimates by con- 
sidering the anisotropic momentum distribution for hard 
partonic degrees of freedom has been reported in Ref . Il9l| . 

At this juncture, we note that the estimates of [l^ 
for the dissociation temperature in the above analyses is 
based on the assumption that the QGP medium is fully 
equilibrated and that it can be treated as an ideal gas 
of quarks and gluons. On the other hand, one knows 
that the latter assumption is particularly far from being 
true. We are also in possession of equations of state for 
finite ternperature QCD which have been evaluated up 
to 0(5^)[2i, [13 and 0{gHn{l/g))^. It is, therefore, 
natural that we should revisit the problem where the 
strongly interacting nature of the quarks and gluons in 
their deconfined state is explicitly and fully incorporated 
in determining this important quantity. Taking into con- 
sideration that the QGP which is produced in heavy ion 
collisions does not possess the full isotropy, but only a 
cylindrical symmetry at most, it is also not out of place 
to study screening phenomenon in an anisotropic QGP 
medium. This sets the motivation for our investigations. 

In the present paper, we shall consider an anisotropic 
QGP (non-equilibrium environment) and employ a re- 
cently proposed quasi-particle model for hot QCD equa- 
tions of state illsi. 
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We argue that all the assumptions taken in Ref.|T9| 
will be equally valid in the present case too. Here, the 
anisotropy is in the momentum distribution functions of 
the hard partonic degrees of freedom same as in Ref . [l^ . 
This situation may be expected in the initial stages of the 
collisions due to the Bjorken expansion and extensi vely 
studied. There is a wealth of literature present |3ll - [3a | 
where the impact of the anisotropy in various observ- 
ables for QGP has been well studied. In most of these 
studies, the isotropic distribution function have been 
considered as the combination of the ideal Bose/Fermi 
distributions [37| . The assumption that the hard degrees 
possess anisotropic momentum distribution functions, is 
only justifiable at the late stages of the collision, so that 
soft degrees of freedom have already get time to ther- 
malize. As correctly stated in Ref.[13, this does not lead 
to a self-consistent non-equilibrium ground state and the 
system has techyonic modes. Nevertheless, as empha- 
sized in Ref.[l^, these issues will be absent in the small 
anisotropy limit. Therefore, in view of this, we only re- 
stricts in the small anisotropy limit. 

Here, we shall consider the more realistic distribution 
functions which are extracted from the improved pertur- 
bative QCD equations of statejl^ [s^]. In this manner, 
we include the interactions in the various observables for 
QGP in the presence of anisotropy. More precisely, we 
address questions such as, (i) the form of the heavy quark 
potential with the realistic equations of state, and (ii) the 
dependence of the dissociation temperature of quarko- 
nia states on the anisotropy and interactions. We shall 
see that the interactions non-trivially modify the static 
heavy quark potential and the dissociation temperatures. 

The paper is organized as follows. In section II, we 
shall review a quasi-particle description of hot QCD EOS 
which will be employed in the subsequent sections to 
study the fate of quarkonia in an anisotropic QGP. In 
section HI, we shall discuss the heavy quark potential 
and melting of heavy quarkonia states. In section IV, we 
shall conclude the present work. 



II. QUASI-PARTICLE DESCRIPTION OF HOT 
QCD 

In this section, we briefly review a quasi-particle de- 
scription of hot QCD equations of state proposed in 
Refs.[29l. Isoj . In these studies, a quasi-particle model 
has been developed to capture all the interaction effects 
present in hot QCD equations of state in terms of non- 
interacting quasi-gluons and quasi-quarks having temper- 
ature dependent effective fugacities. The model came out 
to be quite successful to yield the hot QCD equations of 
state accurately. The model has further been employed 
to study the chromo-electric response functions [301 
transport properties of QGP in RHIC [H, HI]. In this con- 
text, pure SU (3) lattice EOS has also been studied [39j. A 
very strong prediction of the model is regarding the trans- 
port properties of QGP. In particular, the model has been 



employed to study the shear viscosity, 77 of QGP[38l[39|. 
We obtained rj by emplmdng the model in an appropriate 
transport eauation[38ll40(. We showed that both rj and 
its ratio with entropy density, 77/s are very sensitive to 
the interactions and they could be thought of as good di- 
agnostics to distinguish various EOS at RHIC and LHC. 
A very small value for the ratio 77/s was obtained. We 
found that rj/s may violate the universal bound conjec- 
tured from AdS/CFT studies[4l|. 

Here, we consider the quasi-particle description[29l Isoj 
of 0{nh hot QCDIH, [23 and 0{g^ln{l/g) hot QCD 
EOS [ii]. We denote former EOS as EOSl and latter as 
E0S2. EOSl and E0S2 qualitatively matches with the 
more realistic lattice EOS for T > 2 T^. In iacL E0S2 
shows better agreement as compared to EOSl [33] ■ Note 
that EOSl and E0S2 have been evaluated at vanishing 
baryon density. We start with the equilibrium distri- 
bution functions extracted from these two equations of 
state. They posses a very simple form, 
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Zg exp(-/3p) 
1 - 2g exp(-/3p) 

Zg exp(-/3p) 
1 + Zq exp(-/3p) 



(1) 



where g stands for quasi-giuons, and q stands for quasi- 
quarks. Zg is the quasi-gluon effective fugacity and Zg is 
quasi-quark effective fugacity. These distribution func- 
tions are isotropic in nature. At this juncture, we want 
to emphasize the physical meaning of Zg and Zg. These 
fugacities should not be confused with any conservations 
law (number conservation) and have merely been intro- 
duced to encode all the interaction effects at high temper- 
ature QCD. At this juncture, we wish to emphasize that 
the effective fugacities with the similar motivations to 
address com plic ated interactions in BEG (Bose Einstein 
Condensate) [42} and interacting Fermi systems [i^ have 
already been introduced in condensed matter physics. 

Their physical meaning could be understood in terms 
of modified quasi-parton dispersion relations. 



ujf,=p + T^dr Hzg) 
ul^p + T^drHzq), 



(2) 



where Wp denotes the single quasi-parton energy. It is 
clear from Eq. ^ that Zg/g non-trivially contribute to 
the single particle energy. This tells us that the model is 
in the spirit of Landau theory of Fermi liquids. Further, 
Zg/q leads to the quasi-parton number densities Ng and 
Nq as follows. 



g/q 



PolyLog[3, ±; 



■9/9J 



(3) 



where the function PolyLog[3,z] = P"- Both Zg 

and Zg have a very complicated temperature dependence 
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and asymptotically reach to the ideal value unity [30|. 
Here, i^g = 2{N^ - 1) and h'q = AN^Nf, are the num- 
ber of gluonic and quark-antiquark degrees of freedom 
respectively for SU{Nc). The quasi-particle description 
does not change the group velocity of quasi-partons, 
vf,^'^ ~ dpujp^'^ = p (this observation may play an im- 
portant role in setting up the kinetic equation to deter- 
mine the transport parameters for QGP with the quasi- 
particle model). This property allows us to distinguish 
the present model from other models based on an effec- 
tive mass description [i^. lisj . 

This description, which enables us to adapt realistic 
QCD EOS to make predictions for RHIC, leads to very 
very interesting results as far as the bulk and transport 
properties of QGP in RHIC concern. Encouraged by that 
we now seek to employ the understanding to study the 
fate of quarkonia in a anisotropic QGP medium. 



III. HEAVY QUARKONIA IN ANISOTROPIC 
QGP 

It is phenomenologically quite interesting to study 
QGP system which is slightly away from equilibrium. 
It is interesting to inspect quarkonia in a hot QCD en- 
vironment where the hard parton distribution functions 
posses anisotropy in their momentum distribution func- 
tion. This phenomenon is inspected in the recent past in 
Ref.fiot and very interesting physical understanding has 
been obtained. Further, this problem has also been con- 
sidered for an ideal EOS in Ref.j20|. These results deviate 
slightly from [l^ because of the fact that the physics 
setting is different. Thevj20| also assumed non-thermal 
distributions for the soft glouns. For similarities between 
these two analyses, we refer the reader to [l9| . 

Here, we shall primarily follow the work of ^19^ and 
extend in the case of EOSl and E0S2 by combining the 
analysis with the quasi-particle description discussed in 
the previous section. The main motivation here is to 
adapt the realistic hot QCD EOS and investigate the 
fate of quarkonia bound states in this interacting envi- 
ronment. 



A. The real time static potential 



d^p 
(27r)3 



(4) 

where v = {l,p/p), and R indicates that the self-energy 
appears in the inverse of the retarded propagator. The 
hard mode momentum distribution has been assumed to 
be the following form[47j. 



f{p) - foiV^TWW), 



(5) 



where ^ a real parameter known as the anisotropy pa- 
rameter and n is the direction of anisotropy. Without 
the loss of generality, h can be chosen to be along the 
z-direction. The quantity p ■ h = Pz, where Pz is the 
momentum component in the beam direction in RHIC. 
The distribution functions f{p) can be constructed from 
an arbitrary isotropic distribution function, /o by the 
rcscaling of only one direction in momentum space. We 
note that ^ is in [—1, oo], and ^ = corresponds to the 
isotropic case. The values corresponding to the Bjorken 
expansion induced anisotropy are ^ > 0. 

Now, making use of the change of variable p = 
\/p^ + £,p^ , and denoting pz = p cos Op, and following 
the analysis in Ref.fl^, we get the following form of the 
component H™ (which is needed for the real time static 
potential) near the static limit assuming the anisotropy, 
^ to be small (first order in f )[T9l|. 



in kf) 



i-icos(20fc) 



1 + ^ cos(26'fc 



(6) 



where Ok is the angle between k and z. The quantity m u 
is defined in terms of the equilibrium (isotropic) distri- 
bution function as, 



-9' 



d^P dfisoip) 



(27r)3 dp 



(7) 



So far the analysis is quite general and allow any suitable 
arbitrary form for the fiso- Note that often, fiso is taken 
to be a combination of ideal Bose-Einstein and Fermi- 
Dirac distribution functions as 1371, 



Here, we combine the quasi-particle description of hot 
QCD described in the previous section to the formalism 
of heavy quark-potential in anisotropic QGPfl^. The 
aim here is to see how the inclusion of interactions com- 
ing from the realistic equations of state for QGP affect 
the heavy quark-potential and quarkonia melting in an 
anisotropic QGP. 

The starting point to obtain the real time static 
quark-antiquark potential is the Hard thermal loop gluon 
self energy, H'^''^ for a general momentum distribution 
function [33, E^, f{p), 



f^so - 2Ncnbip) + 2Nf{ng{p) + ng{p)), 



(8) 



where nf,{p) 



exp{-Pp)-l) 



ical potential Uq 



and in the zero baryon chcm- 
^ This combina- 



— ii-q ~ exp(-f}p) + l) 

tion leads to the leading order HTL expression (m|j = 
g2(r)T2(7Ve/3-hiV//6)) for the Debye mass in hot QCD. 
Here, Nc denotes the number of colors and Nf the num- 
ber of flavors. 

A more appropriate choice for fiso would be a suitable 
combination of quasi-parton distribution functions [SOj 
(/q and /g displayed in Eq. ([T]).). These are mandated 
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by the fact that QGP is not an ideal gas of gluons and 
quarks rather a strongly interacting system |48l - l53l |. 
We make the following substitutions for fiso, 



J ISO — 



2NcfQ for pure gauge, 
2NJ^ + 2Nff^ for full QCD 



(9) 



Note that the expression for the Debye mass in Eq.© 
with fiso in Eq.(|8|) is same as leading order HTL ex- 
pression obtained from effective kinetic models [s^ - ISTj 
for hot QCD medium with ideal EOS. These results are 
equally applicable in the present case too due the fact the 
the quasi-partons in our quasi-particle model are non- 
interacting. The only consistency requirement is that in 
the asymptotic limit the results should match with those 
for the ideal EOS. This could be achieved with the cho- 
sen expressions for fiso in Eq.®. We shall discuss the 
Debye mass with these choices of fiso in detail in the next 
subsection. 

We now discuss the real and imaginary part of the 
static potential which depends upon the form of mu and 
the anisotropy parameter. In the leading log order the 
real and imaginary parts have been analyzed in Ref. jToj , 
and they obtain the following form, 



Re[V >] 
Im[V >] 



-^[i + o(f)]U^ 

47rr s^J-f 

-g'C,T^H^)[l + 0{-^)] 
3 mor ln[r) 



x(l - -) 



(10) 



where M is the pole mass of heavy quarks |19|. and f ~ 
mor. Here, only the s-wave ground states have consid- 
ered. In this case the term proportional to (1 — 3cos^9k) 
in Eq.® average out to zerofl^ at first order in per- 
turbation theory and corrections will be of the order . 
Similar mathematical forms for the real and imaginary 
part of the static potential are obtained for the present 
case too. The interactions appear via mjy and merely 
renormalizes the effective charges of the quasi-partons, 
which we shall see next. 



1. The Debye mass mo 

Before, we employ the quasi-parton distribution func- 
tions in Eq.(IT]) to obtain the temperature dependence of 
7Ti£), We wish to discuss a very important point here. 
There is an issue of defining the parameter j3 = 1/T' 
which a ppe ars in fiso in terms of the equilibrium temper- 
ature T[19|. The parameter T' will be different from the 
temperature T in the equilibrium. This is due to the fact 
that, in the presence of anisotropy the system is away 
from equilibrium. If we assume the system to be very 
close to the equilibrium (case of very small anisotropy) 
then the parameter T' could be taken to be T. As dis- 
cussed in [l9| . the parameter T' can also be fixed from 



isentropic condition. This is quite tedious here due to 
the non-trivial temperature dependence of the effective 
fugacities, Zg and Zq. We do not consider this particular 
case here and stick to the near equilibrium condition. 

Now, employing /q and /q in Eq.([7]) we obtain, mu in 
the purely gluonic case and full QCD as. 



27-2/':^ ^ QPolyLog[2,Zg\ 



for pure gauge 



N, ^ 6PolyLog [2,Zg] 



Nf^ -l2PolyLog[2,-Zgy ^^^^^^^^^^^^^_ 



(11) 

Here, g{T) is the QCD running coupling constant, Nc — 3 
{SU{3)) and Nf is the number of flavor, the function 

PolyLog[2,z] having form, PolyLog[2,z] = X^feli p-- 
These are the same expression which we have obtained 
from the chromo-electric response functions [soj for the 
interacting QGP. The extra factors appearing due to the 
inclusion of interactions can be attributed to the charge 
renormalization in QGP mediumfs^]. Define the effective 
charges as, 

2 _ 2 ^Pol'yLog[2,Zg] 
— 9 2 



2 _ 2 ~12f o^^/^og[2, -Zq] 

— g 9 



(12) 



The expressions for the Debye mass can be rccasted as. 



Q^T^^ for pure gauge. 



for fuU QCD 



(13) 



Both and acquire the ideal value g^ only asymp- 
totically. Moreover, {Ql,Qq} < g^{T)- In other words, 
these expressions for the Debye masses reach their ideal 
values (leading order HTL expression) only asymptoti- 
cally). To see this, we plot the ratio /i^ = jTi^/m^, 



^ g ^ 6 g 
(1 + ^) 



for pure gauge, 
for fuU QCD 



(14) 



as a function of temperature employing the temperature 
dependence of Zg and Zq from Ref.([30|) in Eq. ([T3|) for 



pure gauge theory and full QCD respectively in Fig. 1. 

The quantity ^.^ scales with T /T^ coming from the scal- 
ing properties of Zg and ^gjs^]. The ratio, /x^j is relevant in 
investigating the dissociation temperatures for quarkonia 
states, which we shall see next. 



B. Quarkonia dissociation 

Here, we choose the following criteria for the quarkonia 
dissociation. The temperature at which the magnitude 
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FIG. 1: Behavior of fid as function of T/Tc- The left panel shows the results for EOSl and right panel for EOS2. In both the 
cases, fid reaches the ideal value unity only asymptotically, and fid < 




FIG. 2: Behavior of Td /T[, as function of T/Tc. The left panel shows the results for EOSl and right panel for EOS2. In both 
the cases, it reaches the ideal value unity only asymptotically, and always > 1. 



of the real part of the static potential (real part gives the 
binding energy) displayed in Ea.([T(Il). equals the magni- 
tude of the imaginary part of the potential (imaginary 
part gives the width of the state) is identified as the dis- 
sociation temperature. 

Employing Eq. ([TO|) . We obtain the following condition, 

--g^^^H X 15 

r 3 m£,r i 

Where r - l/g'^M^. Using T = T and Dcbye masses 
given in Ea. (|ll[) . we obtain the following expression for 



the melting temperature in pure gauge theory, 

= X f i In(^) + ln(i)V \ (1 + i). 

(16) 

On the other hand in full QCD, 

g'^/^M /I 1 W-^\V^ 

^ ^ 31/3(09 + ^Q±)k \2 \9± + !^Q±y J 

x(l + f). (17) 

The corresponding ideal expression for pure gauge the- 
ory reads. 
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ln(-)) ' x(l + |) 
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9' 



(18) 



On the other hand in full QCD, 



anisotropy arc same as that for an ideal EOS in[l9|. This 
is not unexpected because the quasi-parton distribution 
functions employed in this work only differ from the ideal 
Bose/Fermi distributions by purely temperature depen- 
dent suppression factors (effective fugacities, Zg, Zq) oth- 
erwise they are non-interacting. 

IV. CONCLUSIONS 



.4/3^ 



31/3(1 + Ar^/6))^ 
x(l + f). 



2^"^(l + 7V//6)' 



ln(- 



(19) 



Note that Egs. p^ . pT)) posses same ^ dependence as 
in[l9j. There the role played by ^ in modifying the po- 
tential as well as the quarkonia melting temperatures will 
be exactly the same. The interesting observations arc 
regarding the inclusion of the realistic EOS (EOSl and 
E0S2). As it is clear from Eqs.([ini) and ^\ that in- 
teractions non-trivially modify the s-wave ground states 
quarkonia melting temperature. To see the effects of in- 
teractions more closely and quantitatively, we consider 
the ratio Td/T^ in the case of both pure gauge theory 
and full QCD. It is straight forward to see that in the 
weak coupling limit (g << 1), the ratio can be given as, 



rpl 



j-Qs -J 1/3 -f-Qj. pypQ gauge, 

' -1/3 



i+- 



for full QCD 



(20) 



Note that the ratio Td /T^ = ^i^^ . Clearly from the 
Fig.l, the ratio is always > 1. It reaches to unity only 
asymptotically. Therefore, the dissociation temperatures 
with non-ideal EOS will always be higher as compared 
to those for their ideal counter parts. This implies that 
the presence of interactions decreases the effect of Debye 
screening m hot QCD. We have plotted the ratio, Td/T^ 
in Fig. 2. Quantitatively, dissociation temperatures are 
~ 20%-6% higher than the ideal value near T ~ 1.5Tc for 
EOSl. On the other hand, for E0S2, these are - 12%- 
5% higher than the ideal values near T ~ 1.5Tc (sec Fig. 
2). In the case of 2- and 3- flavor QCD, the dissociation 
temperatures are more closer to their ideal counter parts 
as compared to the pure gauge theory. This is true in the 
case of both EOSl and E0S2. In other words, differences 
with the ideal EOS are more pronounced in the case of 
pure gauge theory. 

In the anisotropic case with non-ideal EOS, the dissoci- 
ation temperatures vary from the corresponding isotropic 
case by a factor |. Therefore, in the case of non-ideal 
EOS (EOSl and E0S2) for ^ - 1 there is an increase 
of ~ 10% in the dissociation temperatures as compared 
to the equilibrium case. These conclusions regarding the 



In conclusion, we have studied the form of the static 
inter-quark potential and quarkonia dissociation in an 
anisotropic QGP in the small anisotropy limit employ- 
ing quasi-parton equilibrium distribution functions ob- 
tained from 0{g^) and 0{g^{ln{l/g)) hot QCD equa- 
tions of state. We find that the interactions non-trivially 
modify the real time static potential for a heavy quark- 
antiquark pair. In particular, in the leading order the real 
part of the potential remains intact as for ideal EOS. In 
contrast, imaginary part gets non-trivial modifications. 
Interestingly, the dissociation temperatures are also non- 
trivially modified in the presence of interactions. We 
studied both pure gauge theory and full QCD sectors. 
We show that the interactions coming from the realis- 
tic QGP EOS merely renormalize the partonic charges. 
They enter in the potential as well as in the dissociation 
temperatures through the Debye mass. The above con- 
clusions are found to be valid for both the equations of 
states. These results are valid for small anisotropics. For 
larger anisotropics, numerical simulations are required. 
It is clear from our study that for the non-equilibrium 
states at small momentum space anisotropy one can es- 
timate the relative change caused by the anisotropy at 
leading-logarithmic order employing the realistic QGP 
EOS. 

It would be of interest to extend the present study to 
the case of non-vanishing baryon density and recently 
proposed realistic lattice equations of state with physi- 
cal quarks masses [11], and also the chromo-clcctric 
response functions |30l|. Finally, it would be of interest 
to consider the Bjorken hydrodynamic expansion and 
study the viscous corrections to the response functions. 
These investigations will be done separately in near fu- 
ture. A more interesting study would be to develop a 
two-component model which can handle anisotropics in 
both hard and soft partonic modes and their impact on 
the inter-quark forces in hot QCD medium. 

Acknowledgement: VC is thankful to Rajcev Bhalerao 
for many helpful discussions and Nidhi Gour for careful 
reading of the manuscript. He acknowledges Department 
of Physics, IIT Kanpur, India for the hospitality under 
the TPSC initiative. 



[1] T. Matsui, H. Satz, Phys. Lett. B 178, 416 (1986). 



[2] Nora Brambilla, Antonio Pineda, Joan Soto, Antonio 



7 



Vairo, Rev. Mod. Phys. 77, 1423 (2005). [30 
[3] Agnes Mocsy, P. Petreczky, Euro. Phys. J C 43, 77 

(2005); Phys. Rev. D 73, 074007 (2007); Phys. Rev. Lett. [31 
99, 211602 (2007). 

T. Umeda, Phys. Rev. D 75,017502 (2007). 
J. Adams, et al, Phys. Rev. Lett. 92, 112301 (2004); A. [32 
Adare et al., Phys. Rev. Lett. 98, 232301 (2007). 
V. Agotiya, Vinod Chandra, B. K. Patra, Phys. Rev. C [33 
80, 025210 (2009). 

B. K. Patra, V. Agotiya, Vinod Chandra, Eur. Phys. J. [34 
C 67, 465(2010). 

Saumen Datta, Frithjof Karsch, Peter Petreczky, Ines [35 
Wetzorke, Phys. Rev.D 69, 094507 (2004). 

Nora Brambilla, Jacopo Ghiglieri, Peter Petreczky, An- [36 
tonio Vairo, Phys.Rev. D 78, 014017 (2008). 
W. M. Alberico, A. Beraudo, A. De Pace, A. Mohnari, 
Phys.Rev. D 72, 114011 (2005); Phys. Rev.D 77, 017502 [37 
(2008). 

M. Asakawa, T. Hatsuda, and Y. Nakahara, Nucl. Phys. [38 
A 715, 863 (2003); M. Asakawa and T. Hatsuda, Phys. 
Rev. Lett. 92 (2004) 012001. [39 
T. Umeda, K. Nomura, and H. Matsufuru, Eur. Phys. J. 
C 39, 9 (2005). [40 
H. lida et. al Phys. Rev. D 74 074502 (2006). 
A. Jakovac, P. Petreczky, K. Petrov, and A. Velytsky, [41 
Phys. Rev. D 75, 014506 (2007). 

C. Y. Wong and H. W. Crater, Phys. Rev. D 75 (2007) [42 
034505. 

D. Cabrera and R. Rapp, Phys. Rev. D 76, 114506 [43 

(2007) . 

A. Mocsy and P. Petreczky, Phys. Rev. D 77, 014501 [44 

(2008) . [45 
Nora Brambilla, Antonio Pineda, Joan Soto, Antonio [46 
Vairo Rev.Mod.Phys.77, 1423 (2005) and Refs. therein. 
Y. Burnier, M. Laine, M. Vepsalainen, Phys. Lett. B [47 
678, 86 (2009). 

A. Dumitru, Y. Guo and M. Strickland, [48 

[4:9 



larXi v: 0903 .4703 

M. Laine, O. Philipsen, P. Romatsclike, M. Tassler, 
JHEP 0703, 054 (2007). [50 
Y. Burnier, M. Laine and M. Vepsalainen, JHEP 01, 043 

(2008). [51 

M.A. Escobedo and J. Soto, larXiv: 0804 .06911 

M. Laine, Nucl. Phys. A 820, 25c (2009); [52 

„arXiv:0810. 11121 

F. Dominguez, B. Wu, Nucl. Phys. A 818, 246 (2009); [53 
[arXiv: 081 1.1058' 

Chengxing Zhai, Boris Kastening Phys. Rev. D 52, 7232 [54 
(1995). [55 
Peter Arnold, Chengxing Zhai Phys.Rev. D 51, 1906 
(1995). [56 
K. Kajantie, M. Laine, K. Rummukainen, Y. Schroder 
Phys. Rev. D 67, 105008 (2003). [57 
Vinod Chandra, R. Kumar, V. Ravishankar, Phys. Rev. 
C 76, 054909 (2007). 



Vinod Chandra, A. Ranjan, V. Ravishankar, Euro. Phys. 
J C 40, 109 (2009). 

P. Romatschke, Phys. Rev. C 75, 014901 (2007); P. Ro- 
matschke and M. Strickland, Phys. Rev. D 71, 125008 
(2005). 

B. Schenke and M. Strickland, Phys. Rev. D 76, 025023 

(2007) . 

M. Martinez and M. Strickland, Phys. Rev. Lett. 100, 
102301 (2008); Phys. Rev. C 78, 034917 (2008). 
A. Dumitru, Y. Nara, B. Schenke and M. Strickland, 
Phys. Rev. C 78, 024909 (2008). 

R. Baier and Y. Mehtar-Tani, Phys. Rev. C 78,064906 

(2008) . 

A. Dumitru, Y. Guo and M. Strickland, Phys. Lett. B 
662, 37 (2008); A. Dumitru, Y. Guo, A. Mocsy and 
M. Strickland, 'arXiv : 0901 . 1998 [hep-ph] . 
M. E. Carrington, A. Rebhan, Phys. Rev. D 79, 025018 

(2009) . 

Vinod Chandra, V. Ravishankar, Euro. Phys. J C 59, 
705 (2009). 

Vinod Chandra, V. Ravishankar, Euro. Phys. J C 64, 63 
(2009). 

M. Asakawa, S.A. Bass, B. MuUer, Phys. Rev. Lett. 96, 
252301 (2006); Prog. Theor. Phys. 116, 725 (2007). 
P.Kovtun, D.T.Son, A.O.Starinets Phys. Rev. Lett. 94, 
111601 (2005). 

Mingzhe Li, Haixiang Fu, Zhipeng Zhang and Jican 

Chen, Phys. Rev. A 75, 045602 (2007). 

Ji-sheng Chen, Jia-rong Li, Yan-ping Wang and Xiangjun 

Xia, arXiv: 0712. 0205 (cond-mat). 

A. Peshier et. al, Phys.Lctt. B 337, 235 (1994). 

A. Peshier et. al, Phys. Rev. D 54, 2399 (1996). 

S. Mrowczyhski and M.H. Thoma, Phys. Rev. D 62, 

036011 (2000). 

P. Romatschke, M. Strickland, Phys. Rev.D 68,036004 

(2003); Phys. Rev. D 70, 116006 (2004). 

J. Adamds et al., Nucl. Phys. A 757, 102 (2005). 

G. Boyd et al, Phys. Rev. Lett. 75, 4169 (1995); Nucl. 

Phys. B 469, 419 (1996). 

F. Karsch, E. Laermann, A. Peikert, Phys. Lett. B 478, 
447 (2000). 

M. Cheng, Phys. Rev. D 77, 014511 (2008); A. Bazavov 
et al, arXiv: 090 3 .43791 

Yasumichi Aoki, Zoltan Fodor, Sandor D. Katz, and 
Kalman K. Szaboa, JHEP 01, 089 (2006). 
Szabolcs Borsanyi et. al, laxXiv : 1005 ■ 350"8l 
larXiv: 1007.2580' 

J. P. Blaizot, E. lancu, Phys. Rev. Lett. 70, 3376 (1993). 
J. P. Blaizot, E. lancu, Nucl. Phys. B 417, 608 (1994); 
ibid, 421, 565 (1994); Phys. Rev. Lett. 72, 3317 (1994). 
P.F. Kelly, Q. Liu, C. Lucchesi, C. Manuel, Phys. Rev. 
Lett. 72, 3461 (1994); Phys. Rev. D 50, 4209 (1994). 
Daniel F. Litim, Cristina Manuel, Phys. Rep. 364, 451 
(2002). 



